Luttinger Sum Rule in Slave-Particle Theories 
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The long-standing belief is that the mean-field-like decoupling procedures applied to the slave- 
particle representations of the problems with strong local interaction violate Luttinger sum rule. 
The number of occupied resonant states is small and typically equal to the deviation from the 
sum rule, shedding doubt on the overall results. It is therefore illustrated, on the example of the 
Emery model for the high-T c superconductors, that, through the consistent application of the mean- 
field procedure to the Hamiltonian and the propagators, the sum rule is restored and the resonant 
band conserved. In addition to the resonant band, the electron spectrum contains large number 
of occupied states at the bare site-energy of the ion with strong repulsion. These results can be 
straightforwardly generalized to the other similar problems. 
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The study of strongly correlated electron systems has 
been one of the most active research fields of solid state 
physics for a long time. In the present context, it is ap- 
propriate to mention Friedel/ Anderson and Kondo im- 
purity problems, including the case of Kondo lattice, the 
(extended) Hubbard models and their t-J derivatives and 
the Emery model of the CuC>2 conducting planes for the 
high-T c rare-earth oxides with its t-J limit. General fea- 
ture of these models is the appearance of the resonant 
state, or the resonant band in the translationally invari- 
ant case, on the top of a large fcrmionic background. 

Many methods used in the treatment of strongly cor- 
related systems introduce auxiliary fermions and bosons 
jl], |^, [| [|, |5| . Although the transformation to new par- 
ticles is exact at the outset, it actually becomes useful 
in approximations applied to the coupled fermion-boson 
system. Especially simple situation occurs when only 
the on-site interaction is retained and taken to infin- 
ity, because then only one fermion and one boson field 
are sufficient. These two fields are coupled and various 
decoupling schemes, such as mean-field approximations 
with their Gaussian extensions or the local gauge theo- 
ries where the vertex corrections are ignored, lead to rep- 
resentation of the Green function of the physical electron 
as the product (or convolution) of the auxiliary fermion 
and boson Green functions. This results in (satisfactory) 
appearance of the resonant state or band. However, as 
the rule, the weight of the resonance is small and of the 
order of the error in the Luttinger sum rule associated 
with the calculation of the resonance |j| 0, |) . This sheds 
doubt on the ability of the simple decoupling schemes to 
describe the resonant band, its weight being of the order 
of the error. 

The improvements were attempted by including the 
fermion-boson interactions in high orders with vertex cor- 
rections ignored, considering in particular the non-local 



string fields, with the conclusion that the Luttinger sum 
rule cannot be restored @. Alternatively, the time de- 
pendence of the boson fields was postulated unimportant 
in the calculation of the physical particle propagation [^) . 
This restores the Luttinger sum rule, but results in the 
appearance of the unphysical incoherent background in 
the wide energy range of the spectral density. 

A simple argument is hence given that calculating the 
physical particle propagator consistently with the over- 
all decoupling procedure the Luttinger sum rule is re- 
covered without eliminating the resonant band. The 
large spectral background appears, but localized in en- 
ergy to the vicinity of the non-renormalized deep level. 
To be specific, the reasoning is illustrated on the case 
of the 2d Emery model in the slave-boson representation 
0, [12], [b| [IB], treated in the mean-field approx- 
imation, but the reasoning can be straightforwardly ex- 
tended to the other models in similar decoupling schemes. 

In the usual definition of the Emery model [[l6| for 
the Cu02 plane where Ed and e p are the bare Cu and O 
site energies respectively, to is the bare Cu-0 hopping, 
t' is the 0-0 hopping and Ud is the Cu on-site Coulomb 
repulsion. When Ud ~ oo double occupancy is forbidden 
and the original fermion field c R (spin indices are not 
shown for brevity) which annihilates the d-hole on the 
Cu site is expressed in term of boson b R and fermion f R 
as 
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and c R is given as its conjugate. 

The Hilbert space associated with auxiliary fermions 
and bosons is highly redundant, and the physical sub- 
space is characterized by the no-double occupancy con- 
dition 
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Using Eq.(|I]) for copper sites, combined with the 
fermion fields p R associated with oxygen sites, the Emery 
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Hamiltonian transforms to 

s,R i—x,y i—x 7 y 

Y f R bRPR+8,i + h.c.) + t'(p Rx pR ty + h.c.)], (3) 

s 

where R is a site index and S denotes nearest neighbors. 
Total number N of d and p fermions transforms analo- 
gously to H of Eq.(^). The representations of N and H 
are however not unique, as any operator equal to zero in 
Q R = 1 subspace, such as vf R fR,b R b R , v arbitrary con- 
stant, can be added to n R . In other words, n d R can be 
represented as 

n d R = n f R {b R b R + vb R b R ), (4) 

the expression in the brackets equal to one in Q R = 1 
space. Irrelevant in the exact calculations, the choice of 
v in the brackets becomes important for any approximate 
procedure not confined to the Q R ■ — 1 space. Then the 
only way to keep n R in N and H independent on the 
boson variables, as it is in the exact calculations, is to 
choose 

v=-\ (5) 

in Eq.@). 

Noteworthy, Q R commutes with H and N and there- 
fore they are diagonalized in the Q R = 1 subspace. In 
particular, the physical ground state |G) is the state of 
the lowest energy in that subspace. The mean-field slave- 
boson procedure (MFSB) approximates |G) by the prod- 
uct of the ground state \G(bo)) of the displaced harmonic 
oscillator b' = b — &o and the ground state of the n = 1 + x 
free f- fermions |Go(/)} 

\G) « \G MF ) = |G(&o)> x |G (/)). (6) 

\G(b )) has the property that (b R b R ) = {b R ) 2 = {b R ) 2 = 
Oq. \Gmf) has a component outside the Q R ■ — 1 sub- 
space, which is minimized with the requirement that 
Qr = 1 is satisfied at average 

(Gmf\Qr\G M f) = (n f ) +b 2 = l, (7) 

using Eq.(||). Even when \Gmf) is chosen to satisfy 
Eq.(^), its component outside Q R — 1 subspace is still 
present. This can be easily seen from the calculation of 
(n d ) using Eq.(|), (n d R ) = (n f R )(l + (y + 1)6§). From this 
result, it is obvious that, in the presence of long-range 
order bo ^0, the use of \G_mf) should be combined with 
the choice v — — 1 in Eq.(|J) leading the equality 

«) = (n f R ). (8) 

Using further Eq.(^|), boson field can be averaged out 
from H + \J2Qr of Eqs.® and (|). Eq.(|), i.e. Eq.§) 
is essential in evaluating the local term. The result is the 



MF fermion Hamiltonian which describes free f-fermions 
on the CuOy, lattice with the effective (renormalized) pa- 
rameters of the non-interacting Emery model t = foo^o 
and £f = Ed + A, while t' remains unchanged. On the 
other hand, when the fermion variables are averaged out, 
the Hamiltonian Xb' R b' R of the displaced b' R harmonic os- 
cillator is obtained. Here again, the choice v = — 1 is 
important because it eliminates the b' dependence from 
the local term of H + Q R , i.e. it leaves A to govern 
alone the energy of b' boson. 

Noting that the energy of the b' boson field vanishes 
for |G(&o)), t (i.e. bo) and A are found by minimizing the 
ground state energy Eq of the fermion MF Hamiltonian 
at fixed number of fermions n = 1 + x. 

Time-dependent Gaussian fluctuations around MF 
saddle point are considered in the next step. A quan- 
tity convenient to discuss for such a purpose is the single 
particle propagator between the Cu sites. According to 
Eq.® 

G rr ,{t) - (-i){G\Tf R (r)b R (T)f R ,b R ,\G). (9) 

As in the case of rid in Eq.(Q), the choice of operators en- 
tering the autocorrelation function G RR in Eq. (^|) is not 
unique, because in particular v] R ] R b R b R can be added 
to it without changing G d RR 

G RR {r) = (^)bl{G\Tf R {T)f R b R { T )b R +vf R f R b R b R \G). 

(10) 

The f-fermions and b'-bosons decouple, when the f b' 
coupling is omitted in the MF approximation, i.e. G d 
decouples into the product of fermion and boson propa- 
gators. Noting that to this order the b' boson propagator 
is local 

iG RR ,(T)=b 2 {G (f)\Tf R (T)f R ,\Go(f)h R^R', (11) 
and 

iG d RR (r) = (G (f)\Tf R (T)f R \G (f)} 

(G(b )\Tb R (T)b R \G(b )) + v{nl)b 2 , (12) 

where v is kept explicit in order to follow its trace in the 
following discussion. Note that the operator in Eq.(|lO|) 
reduces to n d of Eq.(^) at r = 0. Consistently with 
the MF choice v = —1, Eq.(||) of the MF procedure, 
this choice should be repeated in Eq.(p"2|) when |G) is re- 
placed by \Gmf) in order to satisfy sum rule (n ) = (n* ) 
of Eq.(||). General structure rj 2 — fj 2 of Eq.(Jl^) is then 
obtained for the correlation function G RR when long- 
time/long-range order (bo ^ 0) is present. This structure 
justifies the time-independent choice of the ^-operator in 
Eq.©. 

Mean-field free-particle propagators appear in Eqs.(|ll|) 
and (|l2]). The Fourier transform (R = R' included) of 
the MF v = — 1 free-fermion propagator is 

Gf ^) = r , (13) 

v ' u3- e(k) + it] sign(e(k) - /i) v ' 
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where e(k) is the dispersion of the free f-fermion band and 
|mj(k)| 2 is the probability of finding it on the Cu-site, 
both quantities obtained by MFSB. MF v = — 1 boson 
propagator (G{b )\Tb' {t)^ \G(b )) in Eq.(|l|) is given by 



D b M = 



A + hi 



(14) 



because, as already mentioned, free dynamics of the bo- 
son b' is local and determined at v — — 1 by the fre- 
quency A (Cartesian gauge of the boson field is meant 
here). Combining Eqs.(|l|), (||) and @), G d RR {oS) reads 



0(/x- £ (k))- 



|m/(k)| 



(2tt) 2 ^ v uj — (e(k) — X) —irj 

+ [G^H- i ^(n / )<5H]6 2 (15) 



The first term in Eq.(p"5j) describes Q the local excitation 
of the dispersionless level of the energy Ed — £/ — A. The 
number of states associated with this level is given by 
the integration of 7r _1 ImG d (w]_ over frequency, and, as it 
can be easily seen using Eq.(|l5|), it is equal to (n^).The 
excitation energy A pc i — e p — Sd of the dispersionless level 
at Ed, large when A is large, is observed in the high-energy 
spectroscopies. 

The second term in Eq.(|l5|) together with Eq.(|l]^ 
forms nearly half-filled resonant band described by 
b^G* (uj, k). Band-dispersion of the resonant band of the 
weight 6g is obviously the one of the MF solution ||] . The 
latter shows that the off-site dynamics is determined by 
the small energy scales: t < t Q (and t' < t in the spirit 
of the Emery model). These scales should be observed in 
low frequency spectroscopies, such as ARPES. 

The integrated contribution of the second term in 
Eq.(|l5|) to n^ 1 lmG d (uj) is It corresponds to 

the contribution of the resonant band b^G^(kuj) to {n d ). 
Together with the contribution of the first term 

and isb^n*) of the third, long-range order term it gives 



(n d ) = (nf)(l + (1 + ^)&q)> which for v = — 1 repro- 
duces the Luttinger sum rule (n d ) — (n?) of Eq.(||). In 
other words, the long-range contribution cancels exactly 
the contribution to (n d ) of the resonant band, without 
removing the resonant band itself. 

Unfortunately, numerous previous calculations for the 
Emery, Hubbard or Kondo problems || fjl pi 9| [l?]] 
amount or are analogous to taking v = — 1 in the MF pro- 
cedure when calculating the energy, but leaving |?| ||, ||] 
v = in Eqs.([lO|- 15) for the propagators, v — leads 
to (n d ) — {nf){l + 6q), i.e. to the breakdown of the 
Luttinger sum rule Eq.(|^), the discrepancy arising from 
the contribution of the resonant band. It is therefore 
important to point out that by simple inclusion of the 
long-range/long-time term in Grr, by choosing v = — 1 
in Eqs.(|lO|) and (|l2"|), solves the problem of the sum rule 
without removing the resonant band. 

Finally, it should be emphasized once again that the 
MF content of the above argument is not essential for the 
restoration of the Luttinger sum rule. The latter relies on 
two assumptions. The first is that the ground state can 
be decoupled into the product of boson and fermion 
states, as in Eq.(H). The second assumption is that the 
long-time order is present in the autocorrelation function 
iGflij, leading to the inclusion of the negative v = — 1 
time-independent term into Eq.(|l2]). The Hamiltonian 
has no explicit role in this argument, except that the 
correlation functions appearing in the calculation of the 
(ground state) energy have to be evaluated consistently 
with G. The present argument applies therefore to the 
entire class of Hamiltonians with large, local interactions 
]l7[, amenable to the slave-particle analysis. 
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